SPECTRAL STATISTICS OF "CELLULAR" BILLIARDS 



BORIS GUTKIN 

Abstract. For a bounded domain f2° C R 2 whose boundary contains a num- 
ber of flat pieces I\, i — 1, . . . / we consider a family of non-symmetric billiards 
f2 constructed by patching several copies of J7 along IYs. It is demonstrated 
that the length spectrum of the periodic orbits in 51 is degenerate with the mul- 
tiplicities determined by a matrix group G. We study the energy spectrum of 
the corresponding quantum billiard problem in Q and show that it can be split 
in a number of uncorrelated subspectra corresponding to a set of irreducible rep- 
resentations a of G. Assuming that the classical dynamics in f2° are chaotic, we 
derive a semiclassical trace formula for each spectral component and show that 
their energy level statistics are the same as in standard Random Matrix ensem- 
bles. Depending on whether a is real, pseudo-real or complex, the spectrum has 
cither Gaussian Orthogonal, Gaussian Symplectic or Gaussian Unitary types of 
statistics, respectively. 



1. Introduction 

According to the Bohigas-Giannoni-Schmit conjecture [Tj the energy spectrum of 
a generic Hamiltonian system with classically chaotic dynamics is distributed in the 
same way as spectra of the standard random matrix ensembles within the same sym- 
metry class. In particular, the spectral statistics of spinless single-particle systems 
with time-reversal invariant classical chaotic dynamics can usually be described by 
the Gaussian Orthogonal Ensemble (GOE), or by the Gaussian Unitary Ensemble 
(GUE), if the time- reversal invariance is broken. For chaotic single-particle systems 
with half-integer spin the spectral statistics are typically the same as in the Gaussian 
Symplectic Ensemble (GSE). There exist, however, a few notable exceptions from 
this rule. It is known, for instance, that the presence of additional symmetries might 
lead to a change of the spectral statistics [21 E]- If the system possesses a discre- 
ate symmetry group H, its spectrum can be split into a number of uncorrelated 
subspectra, where each spectral component corresponds to an irreducible represen- 
tation a of H [HE]. The distribution of the energy levels {E^} within each sector 
depends then on the type of the representation a [HI [7] . For real and pseudo-real 
representations the corresponding spectral statistics are of GOE and of GSE type, 
respectively^ If, on the other hand, a is complex, then the corresponding spectral 
statistics are of GUE type. As a result, even time-reversal invariant systems might 
contain subspectra of GUE type provided H has complex irreducible representa- 
tions [21 E]- Other examples of non-standard spectral statistics are provided by the 

1 I am indebted to C. Joyner, S. Muller and M. Sieber for pointing out to me that pseudo-real 
representations of the symmetry group give rise to GSE type of spectral statistics. 
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Laplacian eigenvalues of certain arithmetic surfaces of constant negative curvature 
[8j [9], as well as by some linear hyperbolic automorphisms of the 2-torus (cat maps) 
[TO] ITT] . Anomalous spectral statistics also appear in some chaotic systems with 
several ergo die components [T2| [T3j. 

From the semiclassical point of view deviations from the spectral universality can 
be always traced to a certain anomaly in the length spectrum of the classical periodic 
orbits (PO). For instance, additional geometrical symmetries imply degeneracies in 
the length spectrum of periodic orbits. Note, however, that even in the absence 
of geometrical symmetries degeneracies in the length spectrum might exist. This 
happens, for instance, in the case of arithmetic surfaces of negative curvature, where 
large multiplicities of the periodic orbits lead to the spectral statistics reminiscent 
of the Poissonian distribution. In the present work we consider another class of non- 
symmetric billiards whose length spectrum of periodic orbits is degenerate. These 
billiards are constructed in the following way. Let fl° be a bounded domain on M 2 
with the boundary dQ° containing I flat pieces (i.e., line segments) r^, i — 
The billiard domain Q is then builded up by taking iV copies of Q° and connecting 
them with the opposite orientations along IYs, see fig. [Th. The resulting domain 
Q is equipped with the flat metric and in some cases with an appropriate choice of 
can be embedded into E jj In what follows we will refer to any such domain f2 
as cellular billiard and to Q° as fundamental cell. It should be noted that the same 
construction can be carried out in any dimension. In particular, starting from a one- 
dimensional fundamental cell one can also construct (quantum) cellular graphs, see 
fig. [lb. It is easy to see that the classical dynamics in Q is intimately connected with 
the classical dynamics in Q°. Specifically, for every periodic orbit in Q there exists 
a corresponding periodic orbit in Q° of the same length. The opposite, however, is 
not always true: a periodic orbit in Q°, in general, gives rise to a number (which can 
be also zero) of periodic orbits in Q. As a result, f2° and Q posses the same length 
spectrum of periodic orbits but with different multiplicities, see fig. [1^. 

Consider now the corresponding quantum billiard problem in Q: 



where the function y? satisfies some boundary conditions at dQ and Aq stands for 
the corresponding Laplacian. Note by passing, that the spectral problem for such 
billiards has previously attracted an attention in connection to the famous question 
of M. Katz [T5] : "Can one hear the shape of a drum?" . It was shown [JH] that starting 
from the same initial domain Q° one can construct in certain cases a pair of non- 
isometric domains Q, Q' such that the spectra of Aq and Aq> coincide (for a similar 
construction of isospectral graphs, see [IZl UH] ) - Here we are rather interested in the 
spectral properties of An for a general cellular billiard Q. Note that, generically, 
Q does not have geometric symmetries. On the other hand, the length spectrum 
of periodic orbits in Q is degenerate and one might suspect that the energy levels 
statistics of A^ exhibit an anomaly. 



2 Note, however, that an embedding of Q in R 2 is not always possible and we actually do not 
require it. 




Any? = Ay?, 
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Figure 1: Construction of cellular (quantum) billiards (a) and (quantum) graphs (b). The 
billiard £1 is obtained by connecting six copies of $7° along its flat sides (coloured lines). 
At each flat part of the billiard boundary d£l we impose either Dirichlet or Neumann 
boundary conditions. 



The main goal of this paper is to give a precise description of the spectral structure 
of Aq, based on the assumption of chaotic dynamics in the fundamental cell Q°. As 
we will show, the situation here is reminiscent of that encountered in systems with 
geometrical symmetries. Namely, Aq can be split in a number of subspectra: 

(1-2) A n = 0A£\ 

a£TZ 

where the sum runs over a subset 1Z of irreducible representations for certain matrix 
group G. Note that in general G is not a symmetry group of the billiard domain, 
but rather a structure group which determines multiplicities of periodic orbits in Q. 
(The exact definition of G and the relevant set TZ of its irreducible representations 
will be provided in the body of the paper.) As in the case of geometrical symmetries, 
the spectral statistics of each sector Aq turns out to be determined by the type of 
the representation a - the statistics are of GUE type, if the representation is complex 
and of GOE, GSE types, if the representation is real or pseudo-real, respectively. 
Furthermore, we show that the spectra of different sectors a are uncorrelated. 

The paper is organised as follows. In Sec. 2 we consider the semiclassical trace 
formula for Aq and express the multiplicities r] 1 of the periodic orbits 7 in Q through 
the traces of certain class of permutation matrices. We then show that r/ 7 combined 
with the phases (resulting from the Dirichlet boundary conditions) can be repre- 
sented as characters of the standard representation for some matrix group G. In 
Sec. 3 we apply the trace formula to obtain spectral correlations of An- We show 
that the resulting spectral statistics are, in general, mixtures of GUE, GOE and 
GSE types of distributions. In Sec. 4 we demonstrate that the spectrum of Aq can 
be split into a number of subspectra corresponding to irreducible representations of 
G and derive a semiclassical trace formula for these subspectra in Sec. 5. Finally, 
the conclusion is presented in Sec. 6. 
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2. Semiclassical trace formula 

Let fl° be a bounded domain on M 2 , with a piecewise smooth boundary dQ° 
containing I flat pieces r i; i = 1, . . . I. We will consider the associated quantum 
Hamiltonian := — (/i 2 /2)Aqo, where h is Planck's constant and Aqo is the 

Laplacian in Q° with the Neumann boundary conditions at dQ°. The spectrum 
{Ej 0) }jLo, Ef ] = h 2 \f/2 of is then defined by the solutions of the following 
eigenvalue problem 

(2.1) -A^f = AfVf, v>feL 2 (n°), ^ 0) |a.o = 0. 

Now, take f2° as the fundamental cell and construct a cellular billiard Q by means of 
the procedure described in the previous section. Note that, in general, the billiard 
boundary dfl contains a number of flat pieces corresponding to unpaired sides Tj of 
the copies of Q°, see fig. [TJ We will study the quantum billiard problem in Q for 
mixed Dirichlet-Neumann boundary conditions at flat pieces of dQ. For the sake of 
concreteness we fix the Neumann boundary conditions at the rest of the boundary. 
Let dQi, i = 1, ... i be the flat pieces of dQ with the Dirichlet boundary conditions 
and let dQ = dQ\ Uf dQi be the remaining part of the boundary, then the eigenvalue 
problem 

— Aifj = Xjifj, (fj G L 2 (Q), 

(2.2) VjldQk = 0, k = l,...£, d n (p j \ m =0, 

defines the Laplace operator and the energy levels {-Ej}°^ ' = h 2 ^j/2 of 
the corresponding quantum Hamiltonian H := —(H 2 /2)Aq. In what follows we will 
consider the spectral density functions for the quantum billiards in Q: 

oo oo 

(2.3) S°\E) =J2 5 ( E ~ ^ d ( E ) = J2 S ( E ~ E ^ 

n=0 n=0 

under the assumption that the classical dynamics in the fundamental cell Q° are 
chaotic. 

The spectral function S 0) (E) = d (0 \E) + d^ c (E) can be split into the smooth 
d~(°\E) and the oscillating part whose semiclassical form is given by the Gutzwiller 
trace formula [TU] : 

(2.4) dM(E) = -Lite Yl Aexp (^(E)) + {°ggK5jr}. 

7ePPO(n°) ^ ' 

Here the sum runs over the set of all prime periodic orbits (PPO) in Q° and S* 7 , Ay 
are the action (including Maslov indices) and the stability factor of 7. In eq. (12.41) 
we singled out the contribution of the prime periodic orbits, as only these orbits 
are relevant for the spectral correlations. The contributions from the periodic orbits 
with a number of repetitions turn out to be suppressed by their large instability 
factors, see e.g., [19] . 

Analogously, one can express the spectral density of states for the billiard Q in 
terms of its periodic orbits. Since each periodic orbit of Q is also a periodic orbit of 
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Q°, the oscillating part of d(E) = d(E) + d osc (E) can be represented as a sum over 
the periodic orbits 7 of the billiard Q°: 



(2.5) d osc (E) = ^Re X 7 A ex p(^7(^)) 



J Contribution from I 
I repetition of p.o. | ? 



where X7 = ( — l)^ 7 ^ with fc 7 being the number of times 7 hits the pieces of the 
boundary with the Dirichlet boundary conditions and ^ 7 G Z + being the multiplicity 
of 7. Note that the forms of (12. 4p and (12.51) are almost identical with a noticeable 
difference of additional multiplicity factors x-y m the last expression. Let us show 
now that these factors can be identified as characters of the standard representation 
for some matrix group. 

Let {T k , k = 1, . . .1} be the set of flat components at the boundary of the domain 
n°. For each such component Tk we define an associated N x N matrix in 
the following way. Let Q®, i = 1, . . . N, be iV copies of Q° which compose the 
billiard Q. Then af^ = 1 for i 7^ j if is connected to through the side Tk, 

of} = 1 (resp. aff = —1) if the boundary component of dQ® belongs to the 
boundary dfl with the Neumann boundary conditions (resp. the Dirichlet boundary 

conditions) and = 0, otherwise. The set of matrices {<J^ k \ k — 1, . . . 1} generates 
then the group G with the multiplication operation given by the standart matrix 
product. In particular, for purely Neumann boundary conditions on <9f2, <j( k ' are just 
permutation matrices and G is isomorphic to a subgroup of the permutation group of 
iV elements. Since G is a matrix group, it admits the standard representation p, such 
that Pij(cr) = (Jij for any a G G. It is straightforward to see that the multiplicity 
factors X7 can be expressed through the charactors of p. Given a periodic orbit 7 in 
fl°, denote T kl T k2 . . . T k -, ki G {1, . . . 1} the (time) ordered sequence of flat pieces of 
dQ° in which the billiard ball flying along 7 hits the boundary, then 

n 

(2.6) x 7 = Tra 7 , ^7 = II ^ G G ' 

i=i 

The above formula can be interpreted in the following way. Any periodic trajectory 
7 in Q° passing through a point x G f2° gives rise to iV trajectories {71, . . .7^} in 
Q. Each trajectory 7^ starts at the point Xj, where Xj G i = 1 . . .N are the 
lifts of x on Q. Since 7 starts and ends at the same point, the set of endpoints of 
71, . . . 7at is, in fact, a permutation of the initial points {x±, . . . xjy}. Specifically, Xj 
is the end point of 7« if (i, j)'s element of <r 7 is nonzero. As a result, by setting all 
nonvanishing elements of 07 to +1 we obtain a permutation matrix of initial and 
final conditions, where the number of units at the diagonal defines the number of 
periodic trajectories among {71, . . .7^}. Furtheremore, having the elements of 
with both negative and positive signs allows us to take in account different boundary 
conditions on dQ which are relevant for the semiclassical trace formula. 

As an example, let us consider two billiards shown in fig. |2] with the Dirichlet 
boundary conditions on the flat parts of the boundary dfl. 
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Example 2.1. For the billiard in fig. [2h the generating matrices are given by 
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The group generated by a a) a b) a c contains 10 elements: G = {oq(t, oio\o G Go}, 
where 

'l 0\ 

and G = Wo,o a ,a b ,a a a b ,a b a a }. 



-o"i 
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Note also that Go is isomorphic to the group of permutations of three elements. 

For the billiard shown in fig. [2b there is an additional connection between first 
and third copies of the fundamental cell and the generators of the group G are given 
by: 
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Figure 2: Pictures of two billiards considered in Example 2.1. Note that in the second 
case the first and the third copy of VLq are connected along the flat pieces of the boundary 
(marked by green colour), as shown by the arrow. 



3. Spectral correlations 
We now proceed with the calculations of the form factor 

/+oo 
R(x/d)e~ 2wiTX dx 
-oo 

for the two-point correlation function of the spectral level density: 
(3.1) R( g ) = ^{d(B + '-)d(E- C -)) -1. 



SPECTRAL STATISTICS OF "CELLULAR" BILLIARDS 



7 



Using the semiclassical expression (12. 5p for the density of states and expanding the 
actions S(E ± |) ~ S(E) ± ^T(E) up to the linear term one obtains 

1 / , . . . a a* „„„ I ^ fa a \ \ s I _ (^7 + ^' 



(3.2) K(r) = -L ^ X 7 XV A^y exp ^I(5 7 - Sy)) 5 (r 



2T H 




where T 7 , Ty stand for periods of 7, 7' G PPO(fi°), and T H = 2nhd is the Heisenberg 
time for f2. Note that the spectral form factor Kq(t) of the quantum billiard in fl° 
can be expressed in a similar way by setting all multiplicity factors x 7 , xy i n ec i- (13.21) 
to one and rescaling T H by the factor 1/N: 



(3.3) K (t) = ^ /^^;,exp Q(5 7 - (7 




In what follows we are going to establish connection between Kq(t) and K{r). 

The semiclassical expressions (I3.2f3.3p can be used in order to calculate K(t), 
Kq(t) perturbatively, as functions of the parameter r. The leading order r 1 con- 
tribution can be obtained using so-called diagonal approximation, where only pairs 
of the same periodic orbits are considered [H]. The next order r 2 term is due to 
the contribution of pairs of periodic orbits with one selfencounter (Sieber-Richter 
pairs) (20]. In the same spirit the higher order terms r n+1 can be obtained from 
the correlations of non-identical trajectories with n selfencounters |21] . Such corre- 
lating orbits can be organised into families according to their topological structure. 
Each family then, includes periodic orbits with close actions which systematically 
contribute into the sum. As a result, the form factor K (r) can be written in the 
perturbative form as: 

00 

(3.4) K (r) = J2r n c n , c n = ^ fcW, 

n=l seSn 

where the last sum runs over the set S n of topologically different structures of 
periodic orbits having n encounters. For generic form of the length spectrum of 
periodic orbits in chaotic systems with time reversal invariance, the contribution 
kn for each structure s 6 S„ has been explicitly calculated in [2T] and shown to 
reproduce RMT result: 



(-2) n 

GOE : ci = 2, c n+1 = i — '—, n>\\ GUE : c x = 1, c n+1 = 0, n > 1; 

n 

GSE : ci = -, c n+ i = — , n > 1. 
2 4n 

To calculate the spectral form factor K(r) for the billiard Q we will assume that 
for long periodic trajectories the multiplicity factors x-y do not correlate with the 
actions 5 7 . It follows then by (13.21) and (13. 3p that 
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(3-5) K(r) = 1 ^(iVr)"<, c' n = £ 

where kn are as in eq. (13.4)) . and 



iV 

n=l seS n 



^i s) = (X 7 X T ') S , 

with the average over all periodic orbits 7, 7' having the same topological structure 
s G «S n of encounters. Since % 7 , X7' are characters of two group elements <r 7 , oy G G 
the above average over periodic trajectories can be substituted with the average 
over the set G^ B ' = {(a, a)} of pairs er, <r G G compatible with the structure s of 
correlating periodic orbits: 

(3.6) v - ) = jm £ 

o-,o-'eG( s ) 

where the normalisation factor |G^| is the number of pairs in 

It follows from eq. (13. 5p that, in order to evaluate K(r) one only need to know 
the coefficients Vn \ Below we show how to calculate T>n^ for a given structure s of 
the correlating periodic orbits. 

3.1. Diagonal approximation. For the diagonal approximation the two trajecto- 
ries 7,7' coincide and we have: G^ = {(g,g) \ g G G}. This yields 

(3-7) Pi = p&( ff ))' x(<r) = Tnx. 

By the group orthogonality theorem (see e.g., |22j) it follows then 

(3-8) 2>i = E n «' 

where n a is the number of times the irreducible representation a enters into p. If 
n a = 1 for each a, then Pi is just the number of irreducible representations contained 
in p. 

3.2. Non-diagonal contribution. The second order term in eq. (I3.5P comes from 
the correlations of periodic orbits shown in fig. [31 These periodic orbits can be 
represented as unions of two (directed) stretches: a U b, a' U b', where a is connected 
with b (resp. a' with b') at the encounter region. Note that in the configuration space 
a and a' are running close to each other. The same holds true for the stretches b 
and b' which have, however, opposite orientations. Schematically, it is convenient to 
denote such correlating periodic orbits as 7 = ab, 7' = ab, where the "bar" symbol 
stands for an "opposite orientation". Such structure of 7, 7' implies that the pairs 
of group elements (07,07/) G G^ in eq. (I3.6p can be represented in the following 
form: 

(3.9) 07 = gh, cry = gh~ x , 
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where g, h correspond to the stretches a and 6, respectively. It is, therefore, necessary 
to calculate the following group average: 

(3.10) ^^EE^w 1 )- 

' ' heG geG 

This quantity can be easily evaluated using the group orthogonality theorem. 



(3-H) ^ = t4E E E P$(jB)P$Hh 2 )p$(j9) 

' ' h,geGi,j,k=la,/3£ll(p) 

h E ^E^h E E 



IGI 14 — ' m a 

1 1 a£K(p) heG 



where the indices in the last two sums run over the set of real TZ T (p) and pseudoreal 
H pT (p) irreducible representations entering p. 
















a 




' - — m_ 


i 


b i 













Figure 3: Sketch of a Sieber-Richter pair in configuration space (left) and the corre- 
sponding diagram (right). The two partner periodic orbits 7 = ab, 7' = ab depicted as 
solid (black) and dashed (red) lines follow each other at the stretch "a" , but after leaving 
the encounter region (shown as a (blue) rectangle on the right figure) move in the opposite 
directions at the stretches "6" and "6", respectively. 




Figure 4: Five diagrams corresponding to topologically different families of correlating 
periodic orbits which contribute to the r 3 term of the form factor. The structures depicted 
at figures (c,d,e) appear only for systems with time reversal invariance. 



To calculate r 3 term of K(t) one needs to take in account several different struc- 
tures of the correlating orbits, which are shown in fig. 0} In general all such structures 
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can be separated into two categories: structures which are relevant both for systems 
with and without time reversal invariance and structures which are relevant only 
when time reversal invariance is present. The first category is composed of two "uni- 
directional" structures shown in figs. HK,b. Here all correlating trajectories have the 
same direction at each of the encounters. The second category is represented by 
three "bi-directional" structures shown in figs. Hb,d,e, where correlating trajectories 
have different directions at least in one of the encounters. As we show bellow, the 
expression for essentially depends on the type of the structure s. Let us con- 
sider for example the "bi-directional" diagram shown in fig. 0b. In that case the 
two correlating orbits have the structure 7 = abed and 7 = adeb, respectively. As 
before, we can use the group orthogonality theorem: 



(3-12) Vf = ^ E E xigfheMge^hr 1 ) 



g,h£Gf,e£G 

1 1 g,h,feG aen(p) a 

jk E 4 E x ia W a \h>) = £ 4+ E 4- 

Gr mi, / — ' mi ^— ' mi 

1 1 aeK(p) a g,h£G aGKr(p) a PeK pl (p) p 



It is straightforward to see that the same result holds for all "bi-directional" diagrams 
of order r 3 , i.e., = P3 = T>^ . On the other hand, for the "uni-directional" 
diagram in fig. one has 7 = abed and 7 = adeb. This leads to 



(3.13) = 4i E E X(9fhe)x(gehf) 

' ' g,heGf,e£G 



4EE xto/^r'k-vv 1 ) = E 

1 1 g,h&Gf,eeG aeK{p) a 

The same result = holds for the diagram on fig. Hb. Note that both 
expressions (I3.12f3.13p have the same form with the notable difference of the range 
of irreducible representations a appearing there. Namely, the sum in eq. H3. 13[) runs 
only over real and pseudoreal irreducible representations entering p while the sum 
in eq. (13. 12[) includes complex irreducible representations, as well. 

Furtheremore, using the same approach it is straightforward to see that, for a 
general diagram of order r n the corresponding coefficients T$ are given by 



(3.i4) vr= V -^r, E -^r+ E 1 4 r . 

aeK(p) a ae-Rr(p) a P&z pr (p) v PJ 

for diagrams with "uni-directional" and "bi-directional" structures, respectively. 
Substituting (13.141) into eq. (13. 5p and taking into account that K (r) = K goe (t) 



SPECTRAL STATISTICS OF "CELLULAR" BILLIARDS 



11 



for a generic Q° with chaotic dynamics, we obtain: 




Note, that since Pij(g) are real matrices, for every complex representation v entering 
p the corresponding complex conjugate representation v enters p, as well. As a result, 



the last sum in (13. 15)) can be cast into the form 

(3,6) E <*■.)' Of) *™ &) • 

where the sum runs over pairs of all complex representations and their conjugate 
counterparts. Eqs. (13.15) 13.16ft suggest the following spectral structure of A$> 

Proposition 3.1. Let p be the standard representation of the structural group G 
(as defined in Sec. 2) and let 

p= M W H 

aev. T (p) 0ell pr (p) uen c ( P ) 

be its decomposition into a number of real, pseudoreal and complex irreducible rep- 
resentations. Then the spectrum of An can be split accordingly: 

i) For each real representation a entering n a times p there exists an associated 
GOE-like subspectrum with the density d Q = dC^j and the number of degenerate 
levels na- 
ii) For each pseudoreal representation f3 entering np times p there exists an as- 
sociated GSE-like subspectrum with the dencity dp = d (j^) and the number of 
degenerate levels 2np. 

Hi) For each pair of complex conjugate representations (u,f?) entering n v times 
p there exists an associated GUE-like subspectrum with the density of levels d v = 
d Ojr-) and the number of degenerate levels 2n u . 

In the next section we analyse the origin of this spectral decomposition in cellular 
billiards. 

4. Spectral decomposition 

Before turning to the general case, let us consider, as an example, the billiards 
shown in fig. [2j By eq. (13.71) the leading order term of the form factor can be 
straightforwardly evaluated giving T>\ — 2, T>\ = 1 for the billiards in fig. |2^i and 
fig. [2b, respectively. This can be understood, as an indication that the spectrum 
of the first billiard is composed of two independent GOE components, while the 
spectrum of the second billiard has a single GOE component. As we show below, 
this is indeed so, since for the billiard in fig. |2k it is actually possible to find a 
projection operator P commuting with A n . To construct such an operator, consider 
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a continuous function ip G C 2 (Q) on Q satisfying the same boundary conditions 
as in (12. 2p . Let {^i,^,^} be the restrictions of ip on = 1,2,3. Regarding 
(■0i> if>2) V'a) as the components of the three-dimensional vector, define the new set 
of functions {ip[, ip' 2 , ^3} 011 i — I5 2, 3: 



(4-1) 0- = J^Py^, 



P 



-1 


1 -1 




-1 1 




1 -1 


ion 


ip' on Q, 


see 


that V" 



f2i,i = 1,2,3 are given by tp"s. It is easy 
function on Q satisfies the same boundary conditions as ip. As a result, the map 
P : ip — > ip' defines the linear operator P which acts on the domain Doul(Aq) of 
Aq. Since any solution of eq. (12. 2 p is mapped by P into another solution of this 
equation we have [P, Aq] = 0. Furthermore, the property p 2 = p implies that P is 
the projection i.e., P 2 = P. 

Turning now to the general case, for each representation a G TZ{p) entering p 
define the following N x N matrix: 

(4-2) P (Q) = ^£* (Q W>)- 

By the group orthogonality theorem these matrices satisfy for 
a, (3 G 7Z(p) and Y2 a en(p) P = ^-NxN- Furthermore, it is straightforward to check 
that the projections (14.21) commute with p(a) for all cr G G: 

[p^,p(a)}=0, aeK(p). 

We can now use p( a ^s in order to construct projection operators P a 's commuting 
with Aq. For a given state ip G H := Dom(An), with the restrictions {^i, . . . ip^} 
on Qj, i = 1, . . . N let ip' G H be the state whose restrictions on Qj, % = 1, . . . iV are 
given by 

N 

(4-3) ^ = E& 

With each p^ we associate the linear operation P a which maps ip into ip'. It follows 
from the definition of P Q and the corresponding properties of p^ that P a P/3 = P a ^a,i3 
for any a, G 7£(p), E aG ^( P ) P " = 1 and 

[P a ,A n ] = 0, aeTl(p). 

Using these projection operators we can now split Aq into the direct sum 

(4-4) A n = A£\ a£> := P q AqP q , 

aS7e(p) 

where each A^ acts on the subspace 7i a = PaH, % = © Qg -R.( p ) Ha- 

Now, let us analyse degeneracies in the spectrum of each Aq . To this end note 
that if a enters n a times into p, the projection p^ can be split further into the sum 
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p(a) _ Y™2iPi \ such, that the subspaces fif' := p\ h, h = C N , i = l,...n a are 
orthogonal to each other and remain invariant under the action of p. Furtheremore, 
in this case there exists a group of unitary matrices u^ 1 commuting with p(g) for all 
g G G which mix different subspaces h^'s inside = © • ■ ■ © hn a , but leave 
every vector v orthogonal to intact: vS°^v = v. By using previous arguments, we 
can lift the matrices to the linear operators Pj )Q , U a acting on the Hilbert 

space "H. From this follows immediately that the Hilbert space l-i a can be split into 
the direct sum 1-L a = 0™^ "H iiQ , "H iiQ = P^W , where [P it(X , An] = 0, Pi, a Pj, a = S it j, 
i = 1, . . . n a and there is a group of unitary operators U a , [U a , An] = which mix 
different H^a and leave states from % a i intact if a' ^ a. In its turn this implies 
that the spectrum of each A[^ is at least n a times degenerate. 

Remark 4.1. For any real representation a, the degeneracy of the corresponding 
spectral component is given (generically) by the number of times a enters p. It fol- 
lows, however, from Proposition ^. 1\ that for complex and pseudoreal representations 
there should be additional double degeneracies in the spectrum. Indeed, for each 
pair of complex conjugate representations (u, v) entering p, the set of eigenvectors of 
Aq is mapped into the set of orthogonal eigenvectors of (and vice versa) by the 
complex conjugation operation. Since all eigenvectors of An can be chosen to be real, 
and A[^ must have the same spectrum. For every pseudoreal representation 
(3, there exists a unitary operator A, such that AA — —1, where A is the complex 
conjugate of A and /3(g) = Afi(g)A~ l for any g G G, see [22]. Combining A with the 
complex conjugation operation Cv = v we obtain the antiunitary operator t = AC 
satisfying t 2 = —1 and commuting with /3(g) for all g G G. In its turn, this induces 
the antiunitary operator T acting on H, such that [Aq ,T] = and T 2 = —1. By 
the last property vectors if) and Tip must be orthogonal to each other for any if) G H, 
which implies the double degeneracy of the spectrum of A^ (Kramers' degeneracy), 
see e.g., ^9\. 



5. Trace formula for subspectra 

By the decomposition (|4.4j) the whole spectrum of Aq can be represented as the 
union of spectra of the operators A^ : 



(5.1) Spec(A n )= (J Spec(A 



(ah 



It is therefore of interest to obtain a semiclassical expression for the spectral density 
of each A^ individually: 

(5.2) S a \E) = --ImTr ( P a 1 ) . 

7r V E + is - A n J 

To this end we can use the same approach, as in the case of systems with geomet- 
ric symmetries [I]. The starting point here is the following representation of the 
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projected Green's function: 



(5.3) G a (E,x,x) 



P, 



a E + ie-A { 



n) = / dy (y\P a \x)G(E,x,y), 



where G(E,x,y) stands for the Green's function in the billiard f2. Let x^ denote 
the mirror image of the point x in the domain Q( k \ k = 1, . . . N with x being equal 
to x^ G n( m \ for some m. Using then the definition (14. 2 p we obtain from eq. (15.31) 



N 



(5.4) G a (E,x,x) = G (£,a;M *«) 



fc=i 

N 



\G\ 



o-GG k=l 



In order to calculate the oscillating part d^l of the spectral density 

(5.5) d {a) (E) = d {a \E) + d^l(E) = --Im / dx G a (E,x,x), 

t in 

we can now use the standard semiclassical representation for the Green's function 
G (see e.g., [H]): 

(5.6) G sc (E jX ^\x^) = ^ ^exp(^(E)\ 

where the sum runs over trajectories in £1 connecting x^ to x^ k \ and Ay, Sj stand 
for their stability factors and actions, respectively. Note that the above expression 
can be also rewritten as a sum over closed trajectories j[x — > x] in the billiard 

(5.7) G sc (E,x^ m \x^) = ^ PmA^)^^p(^(E)\ 

■y[x— yx] 

with a 7 being the permutation matrix (12. 6p corresponding to the trajectory 7. Sub- 
stituting now ( 15. 7p into ( 15. 4p and performing saddle point approximation in eq. (15. 5p 
we obtain for the oscillating part of the spectral density: 

(5 ' 8) d -^ = ^ Re iE E X^Wx^-V^exp (^S 7 (£)) 

Using the group orthogonality theorem we can perform summation over a and finally 

get 

(5.9) d -(£) = ^Rei Yl X (Q) K)Aexp(^ 7 (£) 

The leading order of the mean spectral density d^ (E) can be also obtained from 
eq. (15 .4p by the integration of the imaginary part of the Green's function over the 
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points x^ m ' = 

N 



(5.10) J(«)(£) = -^^ X ( "V)E^>) JdztolmG{E,xt k \xW) 



reG k=l 

<7(zG ^ ' 



JV|G| 



-l/2\ 



where d = Area(f2) /4irh 2 is the leading order (Weyl term) of the mean spectral 
density of Q. 

It is worth mentioning that using semiclassical expression for S a '(E) one can 
straightforwardly establish the type of spectral statistics for each sector a. To this 
end, let us consider the form factor for crossover spectral correlations between two 
sectors a, (3: 

with d a (resp. dp) being the mean density of the energy levels (multiplets) in the 
sector a (resp. 0) given in Proposition 13.11 Assuming, as before, that the averaging 
over ? X can be performed independently of the averaging over periodic orbit 
actions, the problem of calculation K^> (r) reduces to the evaluation of the group 
average: 

(5.11) Idx^+X^i*- 1 )) (X W *(«) + X W *(^))) S 



(<r,5)6G( s ) 



where the sum runs over all pairs (a, a) E having the same structure s E S n . 
By the group orthogonality theorem this average is equal to zero if a ^ (3 implying 
the absence of correlations between different spectral components. It is therefore 
sufficient to consider the case a = /3. As has been explained in Section 3, for real 
and pseudo-real representations the average (15. lip is given by (l/m Q ) n_1 and by 
(— l/m a ) n_1 , respectively whenever s is a structure contributing to the n-th order of 
the form factor. Applying then the same arguments as in the derivation of eq. ( 13 .15)) 
yields: 



(5.12) K^\t) 



5 at pn 2 a K GOE (r) if a is real, 

S at p(2n a ) 2 K GSE (r) if a is pseudo-real, 



where the additional factor 2 for pseudo-real representations accounts for the dou- 
ble degeneracy of the spectrum. In the case of complex representations the average 
( 15. lip is given by \{l/m a ) n ~ l for structures s of the uni-directional type and zero, 
otherwise. This immediately implies that only the diagonal approximation con- 
tributes to the form factor which leads to: 

(5.13) (r) = 5 a> pn 2 a K GVE (r) , if a is complex. 
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Note finally, that summing up the (rescaled) contributions from all sectors a of the 
spectrum gives again the form factor (13.151) . 

6. Conclusion 

To summarise, we have shown that with each cellular billiard Q and prescribed 
boundary conditions on dQ one can associate certain structure group G and its 
standard representation p(G). The characters of p{G) determine multiplicities of 
the periodic trajectories in Q, as well as the phase factors entering the semiclassical 
trace formula for the corresponding quantum billiard problem. The main result is 
that the spectrum of the Laplacian An can be split in a number of uncorrelated 
subspectra in accordance with the structure of p{G). Namely, for each irreducible 
representation a entering n a times into p{G) there exists an associated r^-times 
degenerate subspectrum {E^} whose mean level density is proportional to the 
dimension m a of a. Furtheremore, for billiards with classically chaotic dynamics 
the spectral statistics of {E^} are of the GOE type if a is real, of the GSE type if 
a is pseudo-real and of the GUE type if a is complex. 

It is worth recalling, that the above spectral structure is reminiscent of the spectral 
structure for systems with a group of geometrical symmetries H, where spectrum can 
be split in accordance with all irreducible representation of H. However, one should 
be cautioned to take this analogy too literally, since in the last case, for instance, 
the spectral degeneracies are determined by the dimensions m a of the irreducible 
representations rather than by multiplicities n a (which are not even defined in this 
case) . 

It is natural to inquire about the connection between the geometrical structure of 
Q and the spectral structure of the corresponding quantum billiard. For a generic 
cellular billiard with a large number of connections between its cells it can be ex- 
pected that the matrix group G generated by N x N matrices <Ji,i = I,... I is 
typically maximum possible matrix group G max for a given N. For exclusively Neu- 
mann boundary conditions on dQ, G max is the group of N x iV permutation matri- 
ces. In this case p(G max ) contains precisely two irreducible representations, implying 
that Spec(— An) is composed of two independent subspectra. For the Dirichlet (or 
mixed) boundary conditions on dfl, G max is the group oi N x N matrices having 
the structure of permutation matrices whose elements <7jj = ±1 take all possible 
combinations of positive and negative signs. It is easy to check that in this case 
p(G max ) is an irreducible representation itself and, therefore, no subspectra appear 
in Spec(— An). On the other hand, for some specific structures of Q, and boundary 
conditions on dQ, a richer spectral structure might, in principle, arise. Clearly this 
happens when a cellular billiard posses some geometrical symmetry. In this case 
the standard representation p(G) must contain a non-trivial number of irreducible 
representations of G. 

One can wonder, whether it is possible to have a non-trivial spectral structure 
of Q without having any geometrical symmetry in the system. The answer to this 
question is positive, as can be seen from an obvious example of billiards with the 
Neumann boundary conditions. Any such billiard Q contains as subspectrum the 
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Neumann spectrum of its basic cell 0°. In fact, one can use this property in order to 
obtain non-symmetric cellular billiards Q with an arbitrary number of subspectra. 
Indeed, it is always possible to construct cellular billiards Q' with the Neumann 
boundary conditions having some symmetry H . The spectra of these billiards are 
composed of a number of independent subspectra corresponding to the irreducible 
representations of H . Taking then any such domain Q' as a fundamental cell, one 
can construct a (larger) cellular billiard Q with the Neumann boundary conditions 
which has no geometrical symmetries at all. This billiard, however, will contain a 
number of independent subspectra provided by the Neumann spectrum of Q f . 

The simple arguments above demonstrate that, in principle, it is possible to con- 
struct non-symmetric cellular billiards whose spectrum is composed of several com- 
ponents. It would be of interest to investigate what kind of spectral structure might 
appear in a general case. One interesting question in that regard is whether there 
exist cellular billiards whose spectrum consists of only GUE or GSE components. 
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